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SUMLIARY 



An analysis is made to siraplafy pressure-drop calcu- 
lations for ncnadlabatic and adianatic friction flow of 
air in smooth cylindrical tubes when the density changes 
due to heat transfer and pressure drop are appreciable. 
Solutions of the equation of motion are obtained by the 
use of Reynolds' analorjy between heat transfer and skin 
friction. Charts of the solutions are presented for 
making pressure-drop calculations. A technique of using 
the cbarts to determinfj the position of a normal shock 
in a tube is described. 



INTRODUCTION 



The heat transfer and pressure drop for the flow of 

air in smooth cylindrical tubes may be calculated with- 
out difficulty only when the density changes that ac- 
company'- the changes of temperature and pressure are 
relatively small. In the present paper sn atte.mpt is 
made to analyze flow with large density changes, such 
as occur at ?iigh ?!ach numbers, in order to simplify the 
calculation of such flov/ and to provide a basis for the 
correlation of data. 

The analysis deoends on the validity of Reynolds' 
analogy between heat transfer and skin friction, which 
leads to the equality of the beat-transfer and skin- 
friction coefficients. By rncuTs of this simplification 
the equation of motion can be solved for the nonadiabatic 
case and the stagnation-temp ?ra':ure variation along the 
tube can then be obtrr; led as a function of the Mach 
number; the solution is ^Iven in the form of a chart. 
The adlabatic case, a limiting form of the nonadiabatic 
case, is also- analyzed in some detail and the solution 
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for the presRure variation is given in graphical form. 
Some discussion is included of available experimental 
data (references 1 to 3) order to evaluate the as- 
sumptions that are roade herein and the applicability 
of the solutions. 



•SYMBOLS 

T absolute temperature, '^F + 

X distance in flow direction, feet 

y distance pernendicular to flovi- direction, feet 

u velocity In flow direction, fret per second 

V velocity perpendicular to flov; direction, feot per 
second 

'Btu 

k thermal conductivity, 



seconds foGt^ x °F/foet ■ 
p mass density, slugs per cubic foot 

\i, absolute viscosity, Sii^J 

seconds y feet 

Pr Prandtl number ( aoout 0.7J for laminar flow of air 
at room temperature)* (^Cp/if) 

T_ stagnation temperatTii'6 , F + 

Ctt nondiirensional heat-transfer coefficient 

C„ nond.1mensional i^Vin-f riction coefficient 
F 

D tube diameter, feet 

G m.a33 rate' of flow -oor unit area, _lll?£l 

feet X seconds 

Re Reynolds number (^'^^"^ 

p static pressure, pounds per squrre foot 
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R gas constant, £ po^ -Round s ^^^^^^ ;l716 for air); 

slugs X °F 
also tube radius, feet 

J mechanical equivalent of heat (778 foot-poimds per 
Btu) 

Cp specific heat at constant pressure, 

slugs X °F 

Y ratio of specific heats (taken as l.koh f'or all 
figures) 

M Mach number / ■■ ^ \ 

wall temperature, '^P + 1|.60 
€ nondimensional temperature ratio (Tg/T^^) 
r distance from tube axis, feet 
A tube cross-sectional area, square feet 
Tj tube length, feet 

isentropic m.ass flow per unit area in throat of 



supersonic nozzle, _£iB.£l 



feet^ X seconds 
K constant 
Subscripts: 
c tube center 

0 initial 

1 before shoclf 

2 after shock 

s stagnation 

A single bar over a symbol denotes an average vdth 
respect to cross-sectional area; a double bar denotes 
an average with respect to cross-sectional area with 
pu as weighting factor. 
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RKrNOLDS ' ANALOGY 



Reynolds' analogy between heat transfer and akin 
friction Involves a slinilarity of velocity and tempera- 
ture fields. For stead-^ low-speed flow of air along a 
flat plate, the velocity and teinperature fields are 
given by the slm-olified boundary-layer Oviuations 



and 



6x 6y " PCp 



2 

6u , 6u 11 b -x , , , 

U\ — -I- vt— = (1) 

oy P . ^ ^ 

If the Prandtl nujnbor :Xp/,< is unity and if the 

plcte is' at uniform tent^erature the (T - T^^,) -field 

and the u-field are oirailar, since their different! si 
eq\xations and boandcry condl t3 ons are then identical. 
For laminar flow of air ct ordinary teniporatures the 
value of- the Prandtl number is apt)roxirr.ately For 
turbulent flov, however, the effective values of jj, 
and k are increased bocause of the not turbulent 
Interchange of rnomontun ard hef't, rcsoeotively . The 
Prandtl number for turoulent flow vlll vary somev/hat 
with position because, on tho basis of riiring-len'^th 
theory, the increase in the efiective values'' of the 
thenral conductivity and the viscosity depend on the 
turbulent mixing length and tho velocity iiradient.' 

.'Measurements given in reference }.\ s}iow thst tho 
temperature and velocity prof^'les In the turbulent 
boundary layer along a flat nlate at un-'form tempera- 
ture are nearly the same even though tho boundary 
conditions of the temperature fnd velocity fields were 
not exactly identical during the ir.easareirsnt. The 
measurem.onts indicate that the use of an average Prandtl 
number oi' unity for a turbulent boundary layer should 
not be greatly In error. Tho roasuromonts In refer- 
ence I4. also shov that tlie effect of density gradients on 
the velocity profile may be neglcotod ab moderately 
small temperature differences. 
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Re^molcls' analogy niay be directly extended to lov;- 
speed turbulent flow of air in tubes if the effecta on 
the velocity profile of the pressure gradients due to 
sifin friction and of the density gradients due to heat 
transfer may be di sregai'ded. This assurrptlon is 
justified b"y the good experiiuenfcal agreement of the 
lov'-speed skin-friction and heat-transfer data in refer- 
ence 5 with Reynnlds' analogy. 

The effect of comijressibilT i-,y on the velocity pro- 
file has been measured by Frbssel reference 2), who 
concludes that for subsonic adiataatlc flow of air in 
smooth cylindrical tuues tho velocity oroflle is not 
affected by compressibility. The effect of compres- 
sibility on the tem_perature profile may be seen from 
the following differential equatioix for two-dimiens ional 
tnotion of a comprecsible fluid in a boundary layer; 

Db 2cpj/ PCp 6y^ \^ 2cpjy ■ 

Equation (2) is given by Goldstein in reference 6 and 
attributed' by him to Busemann f reference 7)« 
Pr = 1 and for steady flov/, 

u + V = x- (5) 

Comparison of equations (1) and (J) shov-'S that for com.- 
presslble flow the (T,^ - T,„) -field and the u-field are .. 
similar for identical boundary conditions. For compres- 
sible flow Reynolds' analogy thus postulates the simi- 
larity of the stagnatlon-te.nper ature and velocity fields. 
F'or Incom.pressible flow this sirallarlt^v' reduces to the 
similarity of the (T - T^.) -field and the u-field. 

Reynolds' analogy rcrults in an equality of the 
nondimenpional heat-transfer and skin-friction coeffi- 
cients ^ vmlch is the basis for the subsequent analysis. 
The nondlmenslonal oosfficlents are given by the 
follov>'ing expressions* 
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k 



6y' 




y=0 



pa' 



2 



(5) 



The single bar denotes an average with respect to cross- 
sectional area. Tiie heat- ti'ansf er coefficient is the 



ratio of the heat transfer 
at the wall to the average 



per unit area per unit time 



total energy (relative to 
the energy where - T„,) of the fluid flowing through 

unit crocs-sectional area per unit tirrs. The okln- 
friotion coefficient is the ratio of the shearing 
stress at the vail to the Evorage jnoTn.entirr of the 
fluid flowing through unit cross-sectional area oer 
unit time. The sirrilarity of the (T^,, - Tg) -field 

and the u-fiold gives the following relationship: 







I. J 


y=0 



Pa;.(T,v - T„) 




y~0 



(6) 



From equations (I4.) to (6), 

An experiraental verification of this relationship for 

inconpressible flow is given in reference 5 f-com the data 
of several investigators. 

The value of 0^ will vary with Reynolds nurrher, 

roughness of tube wall, '.;ach rrar.ber, shape of tube 
entrance, and distance fron? tube en1;rance. For odiabstic 
floi.v in smooth tubes beyond the entrance length in which 
the velocity profile *.s still chfjnglng, Frossel in ref- 
erence 2 shows that Cp ic inlependent of i.!ach number 
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•up to the speed of sound and is given by the von Karman- 
Nikuradse relationship 

= -0.8 + 2 log Re v^C^ (7) 



In computing the Reynolds number for the flow in a tube, 
the average air tem.perature may be used for evaluating 
the viscosity. 

Pressure -drop and heat-loss measurements for hot 
air flowing through a water-cooled tube are reported 
in reference 1, For subsonic flov;, it was found that 
Reynolds' analogy Is Independent of Mach number if the 
heat-transfer coefficients are calculated on the basis 
of the stagnation temperature and that the values of the 
skin-friction coefficient ' are in accord with formula f7)' 

Pressure-drop and heat-loss measxarements for burned 
oil-air m.ixtures flov/ing through a smooth water-cooled 
tube at high temperatures and velocities are given in 
reference o. The measurements are in good accord with 
Reynolds' analogy. It is also shown that, under the 
test conditions of reference 0, the skin-friction coef- 
ficients are given by the von F^rraan-Uikuradse relation- 
ship for incompressible flow (equation (7)). 



NONADIABATIC FRICTION FLOW 



The equality of the heat-transfer and skin-friction 
coefficients given by Reynolds' analogy is used to reduce 
the differential equation of motion for compressible flow 
in a smooth tube at uniform wall temperature from three 
variables to two variables. The resulting equation is 
solved by the isocline technique to give a chart for 
making pressure-drop calculations. Particular attention 
is given to averaging the variables in the following 
derivation and solution of the equation of motion: 

The differential equation of miction for uniform- 
pressure at every cross section is 

d (p + pu2) + ^ Cp. pu^ dx = 0 
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The variable pu^ is eliminated by, Introduolng the 
average squared ¥ach number 



='T / dA 



f ri 



^ 2 



0 



pu'^ 
TP 



(3) 



which gives 

d^p + -vpM^) + ^OpYPlv^ dx = 0 
or, in different form, 



I f -i- + 1) dlogp2 + d log ll^ + kc dx = 0 (9) 

The term ^Cp dx Is replaced with a nondiraensional tem- 
perature ratio by equating the differential hoat transfer 
from the wall to the differential heat gained by the 
fluid; that is. 



Pu(T^ - Tg)- CpTTD dx = c^'^ dpuTg 
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D 



d1 



1 - c 



Q 

= li-^ dx (10) 
D 

The double bar denotes an average with respeot to cross- 
sectional area v.-ith pu as 7/eighting factor. introducing 
the last equality of equation (10) in equation (9) gives 



^ ^ + l^ d log p2 + dlog h'^ - d log ( s - 1) = 0 (11) 



2 

The last step in the derivation is to replace log p 

with an expression in ? and F'^. With the aid of the 
relationships 



P = PRT 



Tg = t(i + ^ "^-^ M^) 



T 



2 



2 ' YRT ■ 

p2 is expressed as 

p2 = p2^2 ^—7— N ^^2^ 

= "~2 

The transformation to an expression in e and I.I 
deioends on the considerations that the density and 
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velocity profiles are given by the f- and Ivl^-profiles , 
since 



p = P- 

IiT 



1 + X^ly^ 
2 



and 



These considerations, in conjunction v'ith equations (0) 

— " p 

end (10), show that "c and H" depend only on the €- 
and -profiles, so that equation (12) may be written 

T 

— - r"— r) rp 

P = ^ -7= =A (15) 

v^here the constant K, close to unity, depends only on 
the e- and M^-profiles, or on e and 1,1'^. 

The constant K ho3 been evaluated for the ons- 
seventh power velocity nroflle, v?hich is characteristic 
of fully devolooed tu-.-bulont flow, given by 

u = u, (1 - 
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and for the corresponding slmi'lar stagnation-temperature 
profile, given by 



1 - € 



(1 - 



1/7 



The evaluation of K introduces no particular mathe- 
inatical difficulties but is t_edious. THe values are 

plotted in figuipe 1 against M'^ for several valuer, 
of From figure 1 the variation of K with ¥J^ 

and e is seen to be second order and will be neglected. 
The variation of J\ from %he values of figure 1 for a 
fully developed velocity profile to a value of unity for 
a plane velocity profile is small; therefore K will ■ 
l^ereafter be taken constant nt unity. Taking ¥ con- 
stant is equivalent to assuming one-dimensional flow, 

^ are significan|: averages. 



where e and M 



With the simolif icatlon that 
equation (15) may be written 



F is constant. 



log 



log 



K G^RTv 



+ log € - log 



I?(l + 



Y - 



or, in differential form, 

dlog p2 = dlog e - dlog 



m2 (l 4- I-ZJ. 



(lU) 



Finally, introducing equation (ih) in equation (11) 
gives the_desired differential equation of, motion in 
— p 

e and H , v/hich is 
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Equation (I5) nay be rogardod as a d.lf f erent5 al 
form of the equation of motion when the conditions of 
Reynolds' analogy are fulfilled. Solutions v;era 
obtained by the isocline technique because a closed 
solution was not found. The resulting curves are rhown 
in figure 2. The lower family of curves, for ? < 1, 
corresponds tg heating of the fluid; the \j012er fanily of 
curves, for _c >1, corresponds to cooling" of the_fluid. 
If e >1, e decreases along the tube; and, if c" < 1, 

c increases along the tube. For subson-ic entrance 
I'ach numbers, the ],ocal .?,:ach number' increases vith tube 
length to tho limiting value of unity and, for supersonic 
entrance I.!ach numbers, the local I\ach number decreases 
with tube length to the limiting valiio of unity. In 
oi'der that a solution pass continuously through a Kach 
nurher of unity, the valuo of T would have to increase 
for the upper set of curves, for v-hich the fluid is 
already hotter than the wall, ard would havo to decrease 
for the lower set of curves, for' which tho fluid is 
already colder than tho wall. TnusTuch as such changas 
are contrary to the second law of thennodynaniics , tho 
local f ach number may not -oass continuously through 
unity, 

. Tn figure 2 some er.peri mental data from refer- 
ence 1 are included for comparison wJ th the theoretical 
curves. The tube-wall temoerature was not unlfori?- for 
the experiment; therefore the data haA/e been plotted for 
two separate intorvals of tube length with the tube-- 
wall temperature over esch interval constant at its 
average value. 'i-he agreement between the theoretical 
and experimental results is close. ' Smooth curvos 
instead of eyperlm.ental points are ;!;iv^n in figure 2 
because the experimental points in fij.nres 3 and k of 
reference 1 are connected by smooth curves. 

In order to make more accurate pressure -drop 
calculations for heat exchangers, cliagrDms of f' olotted 
/-•-A 1/2 

against enlarged from figure 2 are given la 

figure 5(a) for subson.tc heating r.nd in fl/T;ure Jfb) for 
subsonic cooling. The cur-'''es of figjrc 5 v.-ere obtained 
bj- ta?-ring the curves of fi^-are 2 as first approximations. 
Special accuracy was attained in rho curves by adjusting 
them so that the integral of the slope along each' curve 
between any two points 
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ciV 



1 n\ 



-.-1/2 



p\l/2 



as given by equation ,■ equaled the difference in 

ordlnatos between the tv/o points on the curve. 

3y equation (15), with Z constant, 

1/2 




(16) 



In order to sirrjplify pressure-drop calculations, lines 
of constant 




1/2 



have been plotted in fij^ure 5(a) and desl (-jnatGd lines of 
constant relative pressure. In rnaking a yressure-drop 
calculation, the value of this expression for the final 
condition is divided by the value for the initial con- 
dition to £lve v/'9o- ^ r^sarao of the procedure follows. 

The value of ?r 



Q for the initial point is found 
from the values of the tube -wall tomperat\ire and 
the initial stagnation teirperature , which either 

ere known or may be Kiessured. Tha value of (M^ ) 
for the initial point is calculatad from equation (IJ) 



for values of G, 



and ID 



o » 



vjh-'ich either are known 



or may be measured. The .values of 



and 



Ik 
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J locate the Initial point in figure 5. The final 

point is located by following the curve of the solution 

,1/2 

through poln1 



' o ' 



U 2 V 



to the final value of e 



This value of e is obtained froia an Integration of equa- 
tion (10), by assuming Cp constant at its average value; 



thus , 



1 + (f^ - l)< 



H D 



(17) 



The values of C^j for fully developed turbulent flow in 

smooth tubes, as given b^~ equation (7), are plotted in 

figure k. For the entrance length of a tube in which the 
velocity profile is changing, the value of Cp' varies 

along the tube and depends on the entrance conditions, 
which determine the initial velocity profile. If experi- 
mental values of for specific entrance conditions 

are not available, the value of Ctt given by figure k 

may be used as the average value for the tube vjith the 
interpretation that this 'value will be slightly Iof. 
Finally, for the initial and final points now located In 
figure 5 the initial and final values of 



M'^/l + ^- 



1/2 



are determined with the help of the lines of constant 
relative pressure, and the pressure ratio p/p,_. is 
calculated by equation (I6). 



As an example, suppose the following values are 
known : 
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S = 0.50 

G = 1.0 



T^, = 760° p abs, 



Re :i 

L _ 

D " 

Pn = 



10 5 
60 

2070 Ib/sq ft 



Equation (I5) gives 



0.5 xH)^ X 1716 y 760 
LIjOJi X (2070)2 



and 



If the average value of Cp for the tube 3 s tal^ron p. 3 
0.00225 from figure 14., the'Value of e computed from 
equation (17) is 



= 1 + (O.5 - l)e 



-Ji(0.00225)(60) 



= 0.709 

Following the curve in figure 5 through the point 
0.525, 0.500 to e = 0.709 locates the point 
corresponding to tube exit. 



The values of 



1/2 



and 



m2 /I + 

\ 2 / 



1/2 



as given by the lines of constant 
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relative pressure are 2.l)j. and 1.66, respectively; 
therefore 



P - 



= 0.775 

The static pressure at the end of the tube is 77.5 per- 
cent of the static pressure nt the tube entrance, 



As previously mentioned, the average value of C 



for the tube is slightly hip:her than the value given by 
figure h; because of the undeveloped flow in the 

entrance lengthy therefore, slightly larger values of £ 
and of pressure drop thtn the foregoing may be antici- 
pated. Curves for constant vslues of 



1/2 



are plotted in figure R ever the rrmge of (_I;fy 
for 'Which an accurate deterinir.atlon is not possible 
from figure 3 • 

Calculations rr.ay be made for tub«s of nonuniform 
v/all tetnperature by porfor.ning the calculation for 
several intervals of the length on tho assumption that 
the tube-wall temperature is unlforiTi' at its average 
value over each interval. The validity of the preceding 
analysis for M > 1 nay not be determined until experi- 
mental data on Reynolds' analogy ore available for 
supersonic flow. 



ADIABATIC PRICTIOII BLO^V 



The preceding analysis for nonadlabatic friction 
flow includes adiabatio friction flow as a special case. 
The equation of. motion lor the adiabatic case is 
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obtained from equation by letting e = 1. The as- 

sumption that e = 1 Is based on the experiinental 
observation that for adlabatlc flow of air the tube 
wall is uniformly at stagnation temperature. Ftxrther, 
equation (3) for the stf gnatlon-ten;perature field Is 
in accord with this observation^ since the particular 
solution 

T = Constant 
s 

corresponds to no heat transfer. 

The reduction of the equation of inotion (equa- 
tion (15)) to the sdiabatlc case proceeds as follows: 

From equation (10), 

(€ - 1) D ' 



Wrxen this expre|.sion 1g introduced in eqiiction (I5), 
together with e = 1, the resulting equation is 



Direct integration gives 



k. f'c^ d21 = 1-/=L- - lr>\+l-tJ:. log. ^ n< 



The integral is retained in the solution because Cv 

may not always be assumed uniform along the tube. The 
variation of C„ will subsequently be discussed. 
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Equation (iG), tlie solution of the equation of 
inotion for adiabatlc friction flow. Is shown in fig- 

urea 6(a) and 6(b) for two subsonic ranges cf M^" and 

In figure 6(c) for a supersonic range of . As 

figures 6(b) and 6(c) show, a particular solution of 
equation (I8) corresponds to a nupersonic value of 

Mq as well as to a subsonic value of Mq • Those 

figures differ only in the lines of constant pressure 
ratio, which will be discussed. 

The solutions are qualitatively like tho se for 

1 < 1 in figure 2. For subsonic values of 1"^'^ the 
local Mach number increases v/ith tube length to a 
limiting value of unity, 'and for supersonic values 

of PIq • the local I.iach number decreases with tube 
length to the limit in;'; value of unit-r. A real flow 
may not oass continuously through lir - 1 because 



is a maximum for this condition. 




. • Before fi.'jure 6 may be used, the variation of 

with Reynolds number, roughness of tube wall, Ivlach 
number, shape of tube entrance, and distance from 

tube entrance should be Icnown. Frossel's data (ref- 
erence 2) show that, for the range of Kach rumiber less 
than unity, the velocity ^profile may not be fully 
developed until L/D = 36; and Kcenan and IJeumann 
(reference 3) found that decreased along the tube 

until l/D was about 56. Beyond th-!^ 3 entrance length, 
Cp may be considered constant. In figvire 7 experi- 
mental data, taken from tho results of refsrencie 5 f'^i* 
flow beyond the entrance length, are compared vjith. the 
theoretical solution of equation (I6), d;, being taken 

constant at the value given by equation (7). The theo- 
retical solution is In good,, agreement with the experi- 
mental data. 



BLGL AHR No. L[|.C16 



19 



The conclusions of references 2 and J for the effect 
of compressibility on Op at supersonic speeds are not 
In agreement. No effect "of Koch number on Op v/as found 
In reference 2 for tubes of L/D < 25, Reference J found 
a decrease In Cp with Mach number for tubes of L/]) ^ ii6 

on the basis of data taken beyond an effective nozzle- 
pipe l/d of about 55 because 'of severe pressure fluctua- 
tions and Dosslble variation of velocity profile In the 
entrance length. Further experiment may be necessary to 
isolate the oompresslblli Ly effect and entrance effect 
on Cp lor supersonic flow. 

Curves of constant pressure ratio p/p^ havfe been 
plotted In figure 6. The oressure ratios have been 
calculated from equation ri6), reduced to 

/.PL\ = ± (19) 

2 

in figures 6(a) and 6(b) the lines of constant 
pressure ratio are based on th e subsonic value of 

¥ ^ and are valid only to M"^ - 1 because the flow may 
0 

not become supersonic for subsonic entrance Ivlach numoers. 
In figure 6(c) the lines of const ant pressure ratio are 

based on the supersonic value of M^"' and are valid for 

the entire range of I^ach mimber. The curves of 
figure 6(0) are also valid in the' case of normal shoc^r 
becaixse the assumptions of equation (19) do not preclude 
shocV. 

As a nui-nerical exam-ole to Illustrate pros sure -drop 

calculations, consider L! = 0.2, Re = ic5, and 
L/D = 60. From figure h, Cxp •- O.OO225. Taking Cp 

constant for the tube as in the example for nonadiabatic 
flov; gives 
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5 ^ = ^5^P 

= ii X 60 X 0.00225 
= 0.5ii 

Prom figure 6(a), p/p^ = 0.755J thus the static pressure 
at the end of the tube is 75*5 percent of the static 
pressure at the entrance of the tube. 



COMrRESSION SHOCKS 



For supersonic entrance VlRoh number-s, snail lengths 
'of rapidly rising pressure may occur within the tube if 
the st.atlc pressure in the reservoir behind the tube 
falls within a certain range. The pressure disturbance 
can be formod by the superposition of a nurtiber of pres- 
sure disturbances orislrating in the dovjnstrea-m flov; 
and travelinp: uostrearn to a point where their velocity 
of propagation is equal to the stream velocity. At this 
point the disturbances would come to rest and form a 
compression shook. Tho position of the shod-: in the tube 
will vary with the static ores sure in the reservoir. 
Figure 6 m.ay be used to determine the positi on , of a 
normal shock with the help of the followin;^ theoretical 
relaticnshiu between thy Liach numbers before and after 



the shock; 



2 _ 2 + - 1) 



= (20) 



2^.1 2 _ ( Y - 1) 



The bars are omitted from the forei^oing equation because 
the shape of the velocity profile in front of and behind 
a pressure disturbance is not known and averages could 
not be evaluated. 

Equation (20) has been 'olottod ' in figure 6. Several 
experimental points from Frossel's results sliow good 
accord with the thocrotlcal curve. The Llach numbers for 
the experimental points have been computed from Prossel's 



nUOU hT(R ¥.0. li+Cl^ 
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data, plotted in terms of G/G^ and p/Pg, by means of 
the following relationship; ^ 



(21) 



where G^ is the isentropic mass flow per unit area in 
the throat of tho suoersonic nozzle in front of the pipe. 
The nondlmensional equation (21) was obtained b^/ lntro-_ 



duclng the following expression for 



in equation (IJ) 




The subscript s refers to the initial stagnation con- 
ditions of the flow. Frossel's pressure rrjeasurements 
were not-made sufficiently close together to determine 
precise shape of the oresGure distribution at the 
For purposes of determining the points of flg- 
a sharp discontinuity was assumed to exist. 



the 
shock. 



ure 8 , 



An exanrale to shovv the method for deternining 
graphically the position of the normal shock is shown 
in figure 6(c). Tf a shock forms at any point in the 
suoer sonic region, the end point of the shock may be 
determined oy moving horizontally to the subsonic Mach 
number g3 van by figure 8. The locus of such end points 



form a shoclr 
in figure 6(c), 



the 



corresponding to a given value of a.^ 

line, which is shown as a broken line 
From the end ooint of the shock, the flow follows 

subsonic solution to the value of h. / Cp dl-j for 

the tube exit; the exit static pressure is then given by 
tho linos of constant pressure ratio. The paths of the 
solution for two positions. of the shock are shown in 



figure 6(c). It 



ma J 



be seen that, c.s the exit presr.ure 



is lowered, the shock moves ilownstream. After the exit 
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pressure is lowered below a particular value, the shock 
v;ill pass out of the tube and pressure equalization will 
occur outside the tube. 

Because the algebraic solutions of the equation of 
motion for adiabatlc flow are given, it is clear that 
the solutions of normal shocks for such flow may also be 
deteraiJned algebraically. Compression shocks in nonadia 
be.tlc friction flow may be handled by the graphical 
technique. Figure 8 raay be used for nonadiabatic flov; 
when the comprossion shock does not occur over too long 
an interval. 
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Figure 1.- Values of K from equation (13) for — -power 
. velocity and stagnation- temperature profiles. 
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- Figure 2.- General isocline solution of equation (15) for 
«'n6nadiabatic friction flow i'n tubes at -uniform, wall - 
/ temperature. ~ -{ Experimental" data from reference 1.) 
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Fig. 3a 



/a 



.6, 



I 

■i . 
I./. 



O. 




i 



tt 



t 



7r, 



77 



± 



i 



i 



z 



7 



!7 



7f 



2: 







Z. 



V 



7 



7 



^2 



cT'' 7 r2 3 ^ B .'y B" 



.9 



;8. 



.7 

.6 
.5 



70 . 



(a) Subsonic heating. ' ' ■ . 

Figure 3.- Isocline solution of equation (15) for air in 
a tube at uni form -wall ' temperature . 
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(b) ^Subsonic cooling. ' - 
- Figure 3.-' Concluded. 
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Figure 4.- Von Karm^n-Nikuradse relationship between 
Reynolds number and skiri-f rlction coefficient for flow 
in sni9ot'h tubes. 
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to' be used with figure 3 for pressure-drpp calculations. 
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Fig.. 6a 
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(a) Low range of subsonic entrance Mach numbers. 

Figure 6.- Plot of, equation (18) for adiabat.ic flow 
. - ' in smooth tubes. 
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-. (b) High range of- subsonic entrance Mach numbers . < 



Figure 6.- Continued. 
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(c) Supersonic entrance Mach numbers. 



Figure 6.- Concluded. ' '' 
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Figure 7.- Comparison ; of equation (18) with experimental, data from- reference 3. 
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Figure 8.- Comparison of equation (20J with experimental data 
- • - ■ from reference. 2. 



